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I. INTRODUCTION 



The electromagnetic transition of the A (1232) resonance to the nucleon A — > A/7* is 
of particular interest |]|||!||| as far as our understanding of the structure of the latter is 
concerned. Historically (e.g. this reaction raised a lot of interest because it allowed 
one to probe the issue of whether the nucleon or its first resonance is "deformed" — the 
reason being that apart from the dominant magnetic dipole (Ml) transition electromagnetic 
selection rules also allow an electric (E2) and a Coulomb (C2) quadruple transition, which 
vanish in simple models of the nucleon with spherical symmetry. Accordingly, the amount 
of deformation can be quantified by the multipole ratios EMR(g 2 ) = E2(q 2 )/Ml(q 2 ) and 
CMR(g 2 ) = C2(q 2 )/Ml(q 2 ), which acquire a four-momentum (squared) dependence in the 
case of virtual photons q 2 ^ 0. 

By the late 1990s there is nearly uniform consensus in the physics community that indeed 
there exist a small quadrupole components in the electromagnetic NA transition (§,§]. In 
the case of real photons one nowadays believes an Re[EMR(0)] « —1/% ... — 4%, but a 
more precise determination of this fundamental property of the nucleon has been surprisingly 
elusive and hotly contested throughout the past decade, both among theorists and among 
experimentalists. It is our hope that the ongoing experiments of electroproduction of 
A(1232) with the resulting better information on g 2 -dependence of the iVA-transition form 
factors, as well as on EMR(g 2 ), CMR(g 2 ) will lead to clear picture of the underlying physics 
and enable us to identify the relevant degrees of freedom for several different regimes of 
momentum transfer. A dramatic change of the physics underlying the electromagnetic NA 
transition is very much expected from perturbative QCD, which predicts for "large" Q 2 = — q 2 
that EMR(Q 2 ) -> +1. At which finite Q 2 the crossover from a negative to a positive EMR 
should happen and whether this point is kinematically accessible at present / future electron 
scattering machines is a further issue of current debate and interest, e.g. M. 

^From a theorist's perspective, treatments of the electromagnetic NA transition may be 
grouped into two categories: 



1. Calculations of the electromagnetic A A/7* vertex with its associated 3 complex form 
factors per se, using different theoretical frameworks that aspire to fundamentally- 
based descriptions of the N — A system. The chiral Bag model(e.g. |J), quark models 
with meson exchange currents (e.g [|K|]), lattice gauge theory (e.g. [l lj) and effective 



chiral lagrangians (e.g. [|T~2| , p~3i| ) constitute such attempts 



2. Direct theoretical treatments |2],|3| of the full scattering processes (e.g. eN — > 
e'Nn, eN — > e'N'j) in the A(1232) resonance region, either based on phenomeno- 
logical Lagrangians supplemented with a method of choice to unitarize the amplitudes 
or dispersion relations. A point of strong contention therein is the issue of separation 
of background vs. resonance contributions. For a recent summary of the status of the 



resulting EMR/CMR extractions we refer to the talk by Workman |14 



In the present work we calculate the A (1232) to N radiative transition A — > A/7* in the 
region of small (i.e. Q 2 < 0.2 GeV 2 ) photon virtuality, utilizing a recently developed effective 
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chiral lagrangian approach [|T|,[nj] that systematically incorporates the spontaneous and the 
explicit breaking of the chiral symmetry of QCD. A small scale e = {p, m w ,8} denoting, 
collectively, small momenta, the pion mass and the Delta-nucleon mass splitting is used 
to establish a systematic power-counting, thus telling us precisely which diagrams/vertices 
have to be included if we want to calculate up to a certain order in e. This approach 
allows for an efficient inclusion of A(1232) degrees of freedom consistent with the underlying 
chiral symmetry of QCD and is referred to as the "Small Scale Expansion" (SSE) ]15|,|16|, 



constituting a phenomenological extension of Heavy Baryon Chiral Perturbation Theory [|T7j 



The formalism can be used to calculate both the vertex as well as full scattering amplitudes- 
here we focus on the former. Clearly, as we are dealing with a low energy effective theory, the 
q 2 dependence of a given physical quantity can be trusted only at low valuesS. Accordingly, 
our final applications will be discussed in this spirit. 

The most general form of the A — > N'y* radiative decay amplitude that complies with 
Lorentz covariance, gauge invariance and parity conservation is described by three form 
factors. We begin our discussion from the widely used formi 



2M N 



M + -^—(Pn -eq^-pN • ? e M ) 



2M 



2M N 



N 



(1) 



Here is the nucleon mass, p^ N A denotes the relativistic four-momentum of the nucleon, 
Delta and g M , are the photon momentum and polarization vectors, respectively. The Delta 
is described via a Rarita-Schwinger spinor w a (pa) with free Lorentz index \i. The Ml, E2 
and C2 multipoles allowed in | — > | electromagnetic transitions can be cast as linear 
combinations of the form factors gi,g2 and g% [ f20|j . 

In this work we study the radiative vertex to 0(e 3 ) in the above mentioned SSE formal- 
ism. This constitutes the first order where pion-nucleon and pion-delta loop graphs enter 
the vertex. As will be shown later, an 0(e 3 ) calculation entails corrections up to (9(1/A^), 
with A x the chiral symmetry breaking scale. One of the main tasks during this calculation 
has been the consistent matching between the results of our perturbative calculation and 
the most general vertex parameterization as given for example in Eq. ([!]). In order to provide 
for a stringent test of our new predictions with experiment we note that a complementary 
calculation of the full pion photoproduction amplitude in the A(1232) region to the same 
order in e 3 is in progress |f2~Tf . 



1 Similar calculations using SSE have been performed for the electromagnetic form factors of the 
nucleon and good agreement with experimental data has been found in the Q 2 < 0.2GeV 2 regime 

I Hill- 

2 Equivalent forms, obtained via use of the equations of motion, can be found in the vast literature 
on this subject. 
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There have been previous analyses of the radiative transition in a similar theoretical 
approach by Butler et al. [12] and by Napsuciale and Lucio [13|. Our work differs from the 
aforementioned references in the following aspects: 



The most crucial difference is that we address the form factors and not just the real- 
photon point. This entails (a) a q 2 dependence in our expressions and (b) the presence 
of an additional form factor (g 3 ), or, equivalently, our calculation yields the CMR(g 2 ), 
in addition to the EMR(g 2 ). 

SSE systematically keeps track of 1/M (i.e., relativistic) corrections to lower order 
couplings. These corrections have not been included in previous analyses |r^JT3|] . We 



also find that our identification of the form factors differs from the one used in |T3[ at 
the real photon point. 



This article is structured as follows. In the next Section we briefly review the essentials 
of the "small scale expansion" and its power counting. In Section III, we discuss the most 
general form of the Pauli-reduced transition amplitude which is consistent with the 0(e 3 ) 
calculation in the SSE. We then proceed, in Section IV, to calculate the loop as well as 
chiral counterterm contributions. The resulting expressions for the form factors are iden- 
tified in Section V, while those for the multipoles in Section VI. Numerical values for the 
EMR(q 2 ), CMR(q 2 ) using presently available (theoretical) input are furnished in Section 
VII. In the same Section we give numerical as well as analytical results for (complex) slopes 
of the three form factors. We summarize and offer our perspective for future efforts in the 
concluding Section. Finally, we devote an Appendix to the discussion of technical matters. 



II. CHIRAL LAGRANGIANS AND THE "SMALL SCALE EXPANSION" 



A. Heavy Baryon ChPT 



QCD being a strongly coupled theory at low (of the order of 1 GeV) energies, renders 
traditional perturbative expansions in the coupling constant inadequate. Chiral perturbation 
theory offers an alternative perturbative expansion, namely one that is realized in terms of 
the external momenta involved in a given physical process. The original strategy was based 
on the notion that at low energies an effective theory of QCD will involve only the nearly 
massless (i.e. the Goldstone bosons: pions etc.) degrees of freedom p2| , p3| . Accordingly, 
chiral perturbation theory has been very successful with respect to applications in the meson 
sector. With the inclusion of baryons, however, the systematic power counting of ChPT 
fails, since baryon masses M^(> 1 GeV) cannot join the set of the expansion parameters 
{external momenta, m^} as they are by no means small and remain finite in the chiral limit. 

A systematic power counting can, nevertheless, be defined through a splitting of the 
nucleonic field degrees of freedom into heavy-light modes and integrating out the former. 
The cost of this procedure is to burden the effective description with additional, higher order 
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contact interactions. Generalizing recent developments in heavy quark effective theories (see, 
e.g., [24 1, ||25|| ) heavy fermion methods were first applied to baryon chiral perturbation theory 
by Jenkins and Manohar The basic premise is the adoption of a non-relativistic mode of 
description, which entails a restriction to four- velocities of the form t> M ~ (1 + \S\, S), \S\ <1. 
On an operational level this means that all momentum dependence in the theory is only 
governed by (non-relativistic) soft momenta k^, defined via 



p M = MoVft + k 



i-i- 



(2) 



where is a typical nucleon relativistic four-momentum and M corresponds to the nucleon 
mass in the chiral limit. The range of validity of the resulting effective theory demands that 
each component k^ « A x , with the chiral symmetry breaking scale being A x « lGeV. This 
approach is commonly referred to as Heavy Baryon Chiral Perturbation Theory (HBChPT). 



In the construction of the effective theory we will follow the systematics laid out in 
The general philosophy is to take the fully relativistic theory as the starting point and then 
perform a systematic non-relativistic reduction. This procedure automatically guarantees 
the proper 1/M corrections to the couplings in higher order Lagrangian terms. Alternatively, 
one could start from a general Lagrangian within the non-relativistic framework, but then 
has to implement the so called "reparameterization invariance" . The latter approach is, for 
example, quite common in the field of heavy quark effective theories 



We now briefly sketch the derivation of the (non-relativistic) chiral Lagrangians for mat- 



ter fields. For details we refer the extensive literature of reviews (e.g. ||16|| , ]T7|, [|29 
We start from the chiral relativistic SU(2) Lagrangian for nucleons 

C N = i) N T N ip N , (3) 
with the relativistic nucleon isospinor field ipjy = ( ^ p J and 

rwff + rf + ... (4) 

being a string of general nucleon-nucleon transition matrices of increasing chiral power 
n 



30[| . For example, to leading order one obtains the well-known structure 



r$ = ip - M + ^ i 75 , (5) 

where denotes the chiral covariant derivative. The parameter g\ corresponds to the axial 
vector coupling constant (in the chiral limit) and the chiral tensor describes the coupling 



of an odd number of pions with nucleon. For more details we refer to |30 
The second step is a redefinition of the relativistic nucleon fields via 
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N = e iMov - x P+ ip. 



N 



H 



AMqv-x 



(6) 



with the velocity-dependent projection operators Pj* 1 = |(1± f)- N is typically called 
the "light" field, whereas H is commonly referred to as the "heavy" field. The relativistic 
Lagrangian Eq.(|3]) then takes the general form 

C N = NA N N + (HB N N + h.c.) - HC N H, (7) 

where the matrices An, Bn,Cn are defined via 

An = Pv^nPv j &n = P v FnPv > Cn = P v ^nP v ■ (8) 



We note that each matrix consists of an infinite string of terms of increasing chiral power n, 
analogous to Eq.([|)!. 

Now one shifts the fields and integrates out the heavy component H . The resulting 
(non-relativistic) effective Lagrangian reads 

C e J f = N(A N + ^BnIoCn'B^N . (9) 



For the case of SU(2) the explicit form of this effective Lagrangian for spin 1/2 nucleons 
has been worked out up to n = 3 by several groups (e.g. [pT,p2,l9[). Generalizations to 
even higher orders are under way [Q. The important point to note is that the inverse of 
matrix Cn is calculated perturbatively, which confines the resulting effective Lagrangian to 
the non-relativistic regime. 



B. The Small Scale Expansion 



HBChPT as described above, has been quite successful for scattering processes off a 
single nucleon near threshold (for reviews see At higher energies, however, the con- 

tribution from nucleon resonances like A(1232) can no longer be parameterized via higher 
order nucleon-nucleon couplings in r^,r^, etc. At some point, i.e. once explicit prop- 
agation of a nucleon- or meson-resonance has to be included, the above described contact 
interaction approach breaks down. So if one is interested in kinematic conditions of such 
dynamic resonance contributions, or in investigations into the low energy structure of nu- 
cleon resonances, it becomes mandatory to include low lying resonances as explicit degrees 



technically speaking, matrices An,Bn start with chiral power n = 1, while Cn begins with 
n = : cl? = 2Mq. This appearance of a large mass term 2Mq > A x in Cn is also the reason why 
the H fields are denoted "heavy" and ultimately get integrated out. 
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of freedom in the effective Lagrangian Eq.(^). In particular this means the inclusion of the 
spin 3/2 nucleon resonance A(1232) in the case of an SU(2) analysis. 



The first efforts towards this direction were performed by Jenkins and Manohar |35 
In the present work, however, we follow a specific generalization of the construction 
method of Ref. [p7]] (as outlined above), which is called the "small scale expansion" (SSE) 
]T3| , P!B|JTP| , PB|J3^7| . The main difference to HBChPT lies in the fact that the chiral power 
counting is modified in a phenomenologically inspired fashion. In the SSE approach one 
expands in the small scale e = {soft momenta, m^, <5o}, where 60 = Ma — M corresponds to 
the Delta-nucleon mass splitting in the chiral limit, whereas in HBChPT one expands in the 
quantity p = {soft momenta, m^jS. The chiral power counting of HBChPT as an expansion 
of all quantities in a power series governed by p n is taken over by SSE as an expansion 
of all quantities in a power series governed by e n . In the following we only briefly discuss 
the construction of the relevant SU(2) Lagrangians with explicit pion, nucleon and Delta 
degrees of freedom and their couplings to arbitrary external fields. For details we refer the 
interested reader to |T3,rB,TU . 



The starting point this time is a set of coupled relativistic SU(2) Lagrangians of rela- 
tivists nucleon and Delta fields iff/ft 

C N = $ N T N tf) N 

£an = r AAr i>N + h.c. , (10) 



where each baryon matrix Tx, X = N, A, AiV is a string of terms of increasing power in e n , 
analogous to Eq.(|). For the spin 3/2, isospin 3/2 Delta field we utilize a Rarita-Schwinger- 



isospurion notation, as laid out in the appendices of Ref. [16|. In the Rarita-Schwinger 
formulation for spin 3/2 fields one encounters the well-known redundancy of degrees of 
freedom, which calls for six projection operators before one is able to isolate the "light" spin 
3/2 component T M from the "heavy" one G^. 

T — p iM v-x p+ p3/2 ,„ / 1 1 \ 



with 



P (33W = 9nu ~ -zlyTtv - ~(Hn v v + •o v rtvfi- (12) 



4 In strict HBChPT 6q counts as a quantity of order p°. This is formally correct but can lead to 
poor convergence properties in the perturbation series. For more details we refer to the discussion 
in |57| regarding the spin-polarizabilities of the nucleon. 

5 For notational simplicity we suppress all explicit isospin indices. 
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The "heavy" spin 3/2 field G M is a five-component object. An explicit representation can 



be found in [fLal but is not needed here. 



Suppressing all spin-isospin indices, the coupled set of relativistic Lagrangians in Eq.(|iO|) 
can then be written as 

C N = NA N N + (HB N N + h.c.) - HC N H 
C A = TA A T + (GB A T + h.c.) - GC A G 
C AN = TAanN + GBanN + HV NA T + HC NA G + h.c. (13) 

Now one integrates out the "heavy" nucleon and "heavy" Delta components H, G and finds 
the resulting (non-relativistic) effective Lagrangian with explicit nucleon and Delta fields: 

C e s f / E = TA A T + NA N N + [TAanN + h.c.}, (14) 



with 



and 



A~a = A A + Io^naIoCn^na + 7o^a7oCa 1 'Ba 
A n = A N + ^B N loC N l B N + iA^qC^Ban 
Aan = Aan + IqVnaIoC^&n + To^AToCA^Aiv (15) 



Cn = C-n — CataC a 1 7oC] va 7o 
Bn = Bn + CnaCa 1 B&n 
T>na = T^na + CnaC^Ba ■ (16) 



The explicit form of the effective Lagrangian Eq. (|14l) in SSE has been worked out to 
n = 2 in Ref. [16]. For n = 1 one finds 



A {1) 

Ana - 
A { P 



-A {1) 

- A/V 
AvA 



D + g A S 

(JkNA 

- (iv ■ D ij - 



u 



s e j 



3/2 



+ 9iS 



u 



9y.v 



(17) 



In addition to the structures already discussed in section |II A| one encounters the Pauli- 
Lubanski spin vector and the 7riVA,7rAA bare coupling parameters g n NA,9i- Further- 
more, the chiral tensors w % ,1$ parameterize the coupling of an odd number of pions to an 
NA, and a AA transition current, whereas D l j corresponds to the chiral covariant derivative 
acting on a spin 3/2 - isospin 3/2 field. With i, j = 1, 2, 3 we denote explicit isospin indices 
and £gy 2 = |5 y — |e ufc r fc corresponds to the isospin 3/2 projector. For more details we refer 
toQ. 



S 



At n = 2 we are only interested in AA^7 vertices, as will be discussed in section IV A . 
The relevant relativistic 0(e 2 ) AN'y contact terms can be parameterized in terms of the low 
energy constants (LECs) b±, as follows 

£atv 7 = Jj^i \9n» + Vil^v] Iplbfi+^N - 2^2^ \9v» + V&l^lA l5fi+D p lp N + h.c. , (18) 

where f?+ = |Tr [/fr ! ] denotes the isovector photon component and y\$ are the so-called 
off-shell coupling constants. 

Performing the transition to the heavy baryon fields for the vertices of interest to our 



calculation one finds |L6| 



M\ = -< f +ia ,S" + ■■■, (19) 



which is independent of the couplings yi,b 6 ,y 6 . In the radiative decay of the A(1232) the 
fact that there is no S-wave multipole allowed precludes the possibility of having an 0(e) 
AA^7 vertex. The parameter b\ therefore denotes the leading 7iVA coupling and carries the 
main strength of the Ml AA^7 transition [r5,IB|,S7|. 



While the desired analysis for n = 2 is quite simple, the situation at n = 3 is more 
complex. In fact, there is no agreement in the literature about number and/or structure of 



the subleading AA^7 vertices (e.g. compare refs. 1 2 , 13LB[ ) . In the next subsection we shall 



explicitly construct the n = 3 AA^7 vertices according to the SSE approach. 



C. Construction of 0(e 3 ) Vertices 



Referring to the SSE master- formula Eqs. flT^ - |T6|) , the relevant Lagrangian contains 
terms which have two separate origins — the dimension six contact interactions contributing 
to A^a an d possible 1/M corrections related to the 0(e 2 ) AN^y vertices. The latter are 
given by 

= N^B^CP-'V^ + ^Bfh^B^T , (20) 

where the superscripts denote the SSE dimension. The pertinent coupling matrices can be 
found in except for the structures denoted as B A j^, T>^\ which only start contributing 



at 0(e 3 ). We therefore construct these matrices, but only for the vertices of interest in our 
calculation. Using Eq. (|T8|) and the orthogonality properties of the light-heavy projection 



operators |T6[ we find (dots signify terms that are irrelevant to the present analysis) 
(2 ) _ ijpi + 2b 6 ) 3/2 ^ + (2) 

NA A liyr ^—^(33) 75Jup V + ••• — U NA 1 + — l /J J 



and (for the b\ part) 
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»(2)6i 
°AN 



4Mn 



,1/2/w 



(i + y)p ( y 2 rv + y/hPllirirf + wP ) 



l/2p,u. 



(21) 

•(1 + 3y)P$'*V + V3yP^(Y + *>) 



(21) 



75/, 



vp 



(22) 



where we have suppressed isospin indices. In a similar manner one obtains the contribution to 
B^' N associated with b&. The complete 1/M corrections at C(e 3 ) which had to be constructed 
via Eg . (|20"D then read 



- (3) fixed 
J A^N-y* 







rpp. 
HP+ i 



4M, 

-^Nf k pP+ S?v%j 2 Di;Tt. 



j 1 



4M 2 



^S 2 S ■ D^Tt 



pp+ v ^3/2 



(23) 



Interestingly, one observes that the off-shell couplings y%, y^ drop out in the final result. 

Finally we have to determine the dimension six jNA counterterms contained in A NA 
as mandated by Eq.([TJp. Once more we have to start from that part of the most general 
relativistic 0(e 3 ) lagrangian that contains all allowed jNA vertices at this order, in which the 



chiral tensor 



|Tr {r* 



Dvi fa/3 



I plays a dominant part. At a first glance the following 

:/3) 9 9 ISJvafii 9 ° l^Jvapi 9 9 I5J va/3i 



structures are possible candidates: g^a a ^ 5 f^ 9^9^15^, 9 m ^lhf^ 9 m 9 vP lsf v + 
e : ^fluBi where the free index fi is to be contracted with the corresponding index of the 
Rarita-Schwinger Delta spinor in the initial state. It is straightforward to verify that only the 
first and third quantities will survive in A^ A and can be shown to give identical contributions 
up to higher order terms. Accordingly we find 



.4 



(3) 
NA 



Pi 

4M 2 



^V^Tr{r* [D v J+ a ]} + 



2M M 



fi OU , AS) | 



(24) 



where Di, E\ are now identified as the 2 7iVA low energy constants of 0(e 3 ). Q This set 
is sufficient for the renormalization of the process A— > N'y*, a more detailed discussion of 
the construction of these counterterms can be found in 



Gathering the above results we determine the relevant SSE 7iVA lagrangian to 0(e 3 ): 



6 The second term of Eq. (|24|) , proportional to E\, appears in the "small scale expansion" due to 
the fact that the Delta- nucleon mass splitting 5q counts as a quantity of 0(e). Therefore one has to 
construct additional relativistic contact terms which are products of lower order interactions and 
the new small scale 5q . For a discussion of this issue see Rcf . [p| . 
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I (3) 



(1) 



T + h.c 



4M 2 



4M 2 



2Mf' Jpf3+ 



(25) 



III. NON-RELATIVISTIC REDUCTION OF THE A — > VERTEX 



Since the leading order Lagrangian counts as 0(e) and, from the general exposition in 
the previous Section follows that each new order in the SSE contributes a power e/M, the 
calculation of the radiative Delta decay to 0(e 3 ) will include all terms up to 0(1/M 2 ). 
Thus, in order to compare our calculation with Eq.(|l]) we must find the most general form 
of the amplitude consistent with an 0(1/M 2 ) calculation. In our case it is convenient to use 
the A(1232) rest frame where a Pauli reduction of the amplitude allows one to identify the 
various multipoles. 

As far as Dirac spinors are concerned the situation with the Pauli reduction is well-known. 
For the Delta field, on the other hand, we think that is worth providing our conventions. 
We start from a typical (e.g. |jj£|| ) Rarita-Schwinger representation of the A(1232) spinor 
with the subsidiary conditions 

7 X = 0, <9X = 0. (26) 



A specific representation consistent with these constraints is the following 

„ f E A + M A \ 1/2 ( 1 \ 



with Ea, Ma and pa the energy, mass and three momentum of the Delta particle, respec- 
tively and S M the four-vector given by 



S • pa 



S + 



S • pa 



M a (£a + M a ) 



Pa 



(28) 



where S denotes the spin 3/2 — > 1/2 transition matrices [B?| . 

Next, we have to expand the various terms entering Eq.(|I|) in powers of 1/M and restrict 
ourselves up to 0(1/ M 2 ) accuracy. Using the analogue of the four-component notation of 
the SSE 



11 



u v (r N ) = P+ u(p N ) 

<A( r A) = P£ M A (Pa) 



(29) 



with Px = Mqv^ + r x , X = N, A, we finally arrive at the following expression 



i M {1) 



|e u 



v(r N ){ (S-e)q lx 



+(S-q)e. 



M, 



N 



+ <D(1/M\ 



gi{q 



6 



5 



+(S ■ q)(v e)q lx 



^(°)-^(°)) +0(1/ M3) 



2M 2 N 



+ (S -q){q- e)q, [o + 0{l/M%)] } < A (0) 



(30) 



which has been evaluated in the rest-frame of the Delta particle. 
Some comments regarding the above equation are appropriate: 

1. To the order we are working the energy u of the outgoing photon four- momentum is 

. . =6 + 0(1/M N ). 



M 2 N + q 2 = 5 



2M N 



2. For q 2 = Eq.(|30"D is not in agreement with the analogous expression of Ref. flTB" 
There, the nucleon mass dependence of the Dirac spinors apparently was not taken 
into account. 



3. Eq.(pO|) has been derived under the assumption that all the form factors are quantities 
of O(e ) at leading order. 

4. It is inconsistent with our subsequent loop results. This can be most easily seen by the 
structure proportional to e • q. According to the 1/M expansion this structure should 
not exist to 0(e 3 ). However, as it will turn out, the loop diagrams of Fig. 1 produce 
a non-trivial q 2 dependence which scales as 1/A?, ~ 1/M 2 . 



We conclude that the popular form of the relativistic isovector nucleon-Delta transition 
current as given in Eq. ([!]), is not ideally suited for microscopic calculations of these form 
factors which rely on 1/M expansions involving a chiral power counting. We find it, therefore, 
advantageousi, in view of the above remarks, to define a new set of form factors, by rescaling 
the previous ones. We propose the following definition 



7 The only alternative would be a g% form factor which behaves as 0{e~ l ) at q 2 = 0, i.e. 53(0) ~ 
However, such a result would violate the basic premise of non-relativistic formalisms like 
HBChPT or SSE that the large scale A x ~ Mjy, 4:7rF v should always appear as A~ n , n > in order 
to ensure a consistent power counting. 
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G 1 (q 2 ) = 9l (q 2 ), G 2 (q 2 ) = g 2 (q 2 ), G 3 (q 2 ) = -^9^) 

M N 



(31) 



which leads to the relativistic A — > N'y* transition matrix element 



Gi(q 2 ) tl , , , G 2 (q 2 ) , 
, G 3 (q 2 ) 



N 



AM N S 



(q ■ e q^ - q\) 



(32) 



Eq.(p2|) will serve as the reference point for future comparisons with other calculations. 
The advantage offered by this parameterization can be seen best when performing a 1/M 
expansion of Eq. (|32|) , once more restricted to O {1/M 2 ): 



iM NR 

tJvl A->N^ 



|e u v (r N 




2M 2 N 



+ 



5Gn 



4M 2 N 



+ 



AM 2 N 5 



G 3 (q 2 ) 



G 2 (q 2 ) 
4M 2 N 



+ 0{l/M%) 



■ig + W) 



<a(0) ■ 



0{l/M%) 



(33) 



The above formula serves as the basic connection between our (non-relativistic) results 
calculated via SSE and the desired (relativistic) transition current Eq. (|32|) . One can easily 
see that to 0(e 3 ) one is sensitive to the leading and subleading behavior of the form factor 
Gi{q 2 ), whereas the form factors G 2 (q 2 ), G 3 (q 2 ) start out at 1/M 2 and therefore only allow 
us to determine them to leading order. One can also check that Eq . (|33|) is gauge-invariant 
to £>(1/M 2 ) and allows for a non-vanishing g 2 -dependent contribution in each of the four 
independent structures, consistent with the SSE calculation to 0(e 3 ). We now move to the 
details of the calculation. 



IV. THE CALCULATION 



A. Born Contributions 



In the present subsection we discuss the Born contribution to the amplitude. The relevant 
expressions resulting from each ANj* tree level vertex have the following explicit forms! 



8 To the order we are working we can identify the parameters Mq,5o with the physical values 
M N ,S = M A - M N . 
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■ Kyfbe rc 



ge u v (r N ) 



- Jfe u v (r N ) 



|e u v (r N ) 



2 e u v (r N ) 



ge u v (r N ) 



q 6l 

5 '^2M 



A/ 



D x 5_ 



2M N 
<a(0) 



<a(0) 



A? 



A? 



c M . c -h 



A? 



AM 



N 



<-> ■ <7 e u 5- + Co <-> ■ q q u s- 



'2M 



N 



N 



<a(0) 
<a(0) 
<a(0), 



(34) 



where the superscripts tabulate the relevant vertices, in relation to the LECs, while "rc" 
stands for "relativistic correction" . 



We mention that, the second relativistic correction Lagrangian term in Eq.(p3|) gives 
zero contribution to our form factors. To see this, note that the derivatives acting on the 
Delta field bring down factors (S ■ r&) and (r A • T) where r A is the A(1232) soft momentum. 
In the A(1232) rest frame, where our calculation is performed one finds r A = 5iA Both 
(S ■ Ta) and (r& • T) therefore vanish due to the light Delta constraints. We do anticipate, 
on the other hand, that this operator will contribute to the resonant Born diagram in pion 
photo-production where the intermediate A(1232) will also contain off-shell components of 
the Delta field. 

The overall Born contribution, therefore, is 



i KA Born 



-e u v lr N 



-61 {2E l - £>i)<T 



2M 



N 



AM 2 N 



+S ■ g e M 



61 



2M 



E x 5 (h + 26 6 ) <f 



N 



2M 2 N 



AMI 



+v ■ e S ■ q q^ 



D 1 bt + 2b 6 



AM% 



AMI 



u 



v.A 



(0). 



(35) 



In the next subsection we will turn our attention to contributions from loop corrections 
to the transition vertex. 



B. Loop contributions 



To 0(e 3 ) one can draw 10 independent loop-diagrams for the A — > A/7 transition (see 
Fig. [T]). However, due to the constraints v ■ u l v A = S ■ A = 7vu* A „ = satisfied by the 
on-shell spin 3/2 isospin 3/2 spinor u l v A/1 only diagrams li, lj in Fig. 1 survive. Their 
contribution can be formally written as 
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iM 



l 2 2eg nNA _ 
3 Fl Ul 



{r N ){{S ■ e)q fl [g A F N + xgiF&] 
+ (S ■ q)e fl [g A G N + X9iG A ] 
+ ( s ■ ?)^( e o \9aJn + XQiJa] 

+{q ■ e) [g A H N + X9iH A }) }<a(0) • 



(36) 



Here the N and A labels on the functions Fj, Gj, ... denote the intermediate baryon prop- 
agator in diagrams (i) and (j) of Fig. 1, respectively, while the isospin factor x takes the 
value — |. The loop functions are defined via 



F N (t) 
FA(t) 
G N (t) 
G A (t) 

Mt) 
Mt) 

H N (t) 
H A (t) 



f 1 , , A 301 (5x,M? 

2 dx (x - 1 — - 

Jo d 



2 / dx 
Jo 

„1 | 30 J 

2 1 dx x 
Jo 



X 



A 301 (-5x, M t 2 ) + 2x A 301 (-5x, Ml 



d 



A 301 (5x,M? 
d 



-2 j dx 
Jo 



x — - 1 — — + 2(1 - x) 

d 

310/ 



d(d- 1) 
A 301 (-5x,M?) 



d(d - 1) 



o 



dx x(l-x)A 3W (Sx,M t 2 ) 



i d — 3 

dx x(l - x)- A 3W (-Sx, Ml) 

o d — 1 

1 dx x(l - x)(l - 2x)A 300 ( y 5x, Ml) 

(i 

1 dx x(l - x)(2x - l)^^-A 300 (-6x, Ml) , 
o d — 1 



(37) 



where the t = q 2 dependence arises implicitly through M t 2 = m 2 + tx(x — 1). Explicit 
representations for the A^ k (fl, Ml) are given in the Appendix. The 1/M 2 sensitivity of our 
calculation is saturated by the two inverse powers of the quantity AttF^ ~ M which is of the 
order of the chiral symmetry breaking scale A x (= {4irF n , M^}). 

We close the present subsection with a brief discussion concerning renormalization is- 
sues. It is straightforward, using Eq.(|37D and Eq. (|A5|) , to see that Eq. (|36D for the loop 
contributions to the decay amplitude of the A(1232) is a complex quantity and contains 
several divergent pieces. The lower-order couplings bi, b G only possess a finite part, whereas 
the infinities encountered in the loop amplitudes are absorbed by the infinite parts of the 
0(e 3 ) countertermsi E x , Di 



9 The interaction terms proportional to these two LECs are what one calls counterterms in the 
conventional field theoretical language. Nevertheless, as is usual in ChPT, we will sometimes 
employ the same name for the b%, be couplings. 
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D 1 = D{(\) + 16tc 2 (3 d1 L 
E 1 = E[{\) + 16tt 2 (3 e1 L . 



(38) 



with 19 



L 



A 



d-i r 



16tt 2 



_l_ + i( 7£ _ 1 _ 111 4,) 



(39) 



The two pertinent beta-functions Pdi, Pei can be found from Eq.(|37|) and the expressions 
in Appendix A. It is straightforward to obtain the results 



fa 



g n NAM%f . . 

9-kNaMn 
6n 2 F 2 (gA-20/9gi). 



(40) 



At the moment, however, we have no information about the numerical size of their finite 
parts E[(\),Dl(\) at the chosen renormalization scale A, as well as about the magnitude 
of the coupling 6 6 . We note that this information will be available in the near future as 
several calculations regarding different scattering processes in the delta resonance region 
using SSE are underway, which will allow for a systematic extraction of the relevant higher 
order iVA-couplings. In section VII we will nevertheless give some numerical estimates for 
linear combinations of these couplings based on calculations for iVA-transition multipoles 
performed in theoretical frameworks outside SSE. 



C. A note on gauge invariance 



Gauge invariance requires that, upon setting e M — > g M , the amplitude must vanish. The 
general amplitudes Eq.(^) and Eq.(|3"3"|) satisfy this requirement explicitly (recall that eo — > 
qo = 5). In the case of our loop calculation this requires 

Fi(t) + Giit) + 8Ji{t) + tHi(t) = , (41) 



separately for i = N and % = A. This provides for a stringent test of our calculation. From 
the explicit formulae Eq.(p7[) and Eq.([A5|) one can show that the constraints of Eq. (|4"1|) , 
are indeed satisfied. This is most easily seen by partial integration of the ^ 2 part of 

(5Ji + tHi). 
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V. THE A -> iV 7 * FORM FACTORS TO C(e 3 ) 

Combining the loop result, Eq. (|36|) , with the Born contributions, Eq.(^), we obtain the 
total 0(e 3 ) amplitude. This is to be compared with the general form of the amplitude, 
Eq. (p3]) . The identification of the three form factors is straightforward. At this point we 
have no information about the magnitude of the counterterms bi,b^, Dx, E,. However, the 
important observation is that the unknown counterterms only affect the overall normal- 
ization of the 3 form factors Gi(q 2 ) and leave the q 2 -dependence unaffected. We therefore 
separate the real photon point q 2 = from each form factor and write 

G 1 (q 2 ) = 0,(0) + ~ gi (q 2 ) 
G 2 (q 2 ) = G 2 (0) + ~g 2 (q 2 ) 
G 3 (q 2 ) = G 3 (0) + h(q 2 ) , (42) 

with <7i(0) = 0, % — 1,2,3. It is then straightforward to arrive at the results of Table I. 
Comparing the relativistic form of the coupling structure proportional to b\, Eq.([L8|), with 
our general relativistic amplitude, Eq.(|32|), one expects that only the Gi form factor receives 
contributions from this coupling, which also holds in the non-relativistic SSE formalism. 
Analogously, the coupling 6 6 in Table | is shown only to contribute to the form factor G 2 , 
again as expected from its relativistic analogue in Eq. (|T8|) . Furthermore, the results of 
Table [0 they should be, gauge- independent. This happens automatically for the 

chiral counterterm contributions. For the loop contributions, however, gauge invariance can 
be explicitly demonstrated via the use of Eq. fl4T|) . 

The results of Table | can now be compared with the ones in Ref. [|X3j] where only the 
real-photon case is considered (no G3). We find a difference from our result originating in 
Eq. (|30D . In particular, there is no contribution to G2 from the b\ counterterm. This happens 
because, in ||13|| , the coefficient of G\ in the (S ■ q)(e ■ T3) term is 2Mn rather than 2Mm + 5. 

In section |V1I| we derive an estimate for the unknown couplings. We will then be in 
position to describe both the q 2 dependence and the overall normalization of the three form 
factors. Once again we stress that chiral symmetry demands that these form factors are 
complex quantities, which is due to the on-shell ttN intermediate state of loop diagram (i) 
in Fig.l. 

VI. CONNECTION TO A -» Nj* TRANSITION MULTIPOLES 

In the relativistic case the identification of the transition multipoles Ml(g 2 ),E2(g 2 ),C2(g 2 ) 
in terms of the "Dirac"-type form factors Gi(q 2 ),G2(q 2 ),G^(q 2 ) is obtained |2(J by reading 
off the spin-space tensor properties of the Pauli- reduced most general amplitude, Eq. (53), 
in the q • e = gauge. The explicit expressions read0 



10 At the photon point (q 2 = 0) one easily verifies that they agree with the corresponding expres- 
sions given in 0. We are grateful to Rick Davidson for valuable support in this point. 
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Ml(g 2 ) = c A {Gi(g 2 )[(3M A + M N )(M A + M N ) - q 2 



E2(q 2 ) 
C2(q 2 ) 

with 



-G2(q 2 )^(Ml -Ml- q 2 ) - G 3 (q 2 )q 2 ^-\ 



2M 



N 



-CA {G 1 (g 2 )(Mi -M 2 N + q 2 ) - G 2 (q 2 )^(M 2 A - M 2 N - q 2 ) - G 3 (q 2 )q 2 ^) 

-c A 2M A \q\ {^(g 2 ) - G 2 (q 2 )^- - G 3 (g 2 )^} , 



(43) 



CA 



e\q\ 



2M A 



12M N M A J2M N (E N + M N ) V m a - M 



(44) 



N 



We now need to match these general expressions to the 0(1/ M 2 ) accuracy of our results for 
the Gi,G 2 ,G 3 form factors obtained via SSE in the previous section. We may proceed in 
two ways: 

1. Use the same reasoning as in the relativistic case, i.e., obtain the multipoles from the 
spin tensor structure of the appropriate Pauli-reduced expression, Eq. (|33|) . 

2. Expand the relativistic formulae for the multipoles to 0(1/ M 2 ). 



In both cases we obtain 
e^ 2 — q 2 



Ml(q 2 ) 
E2(q 2 ) 
C2(q 2 ) 



6VS 
e\/5 2 — q' 2 

e(5 2 - q 2 ) 
6VS 



2Gi( 



M N 
5 

'2Mj 



() G 2 (q 2 ) 



AM 2 N 



A5M 2 N 



Gi(Q) 



AM' 



6 —G 2 (q 2 ) + -1 



N 



G 3 (q 2 ) + 0(1/M%) 
ASM 2 



-G 3 (q 2 )+0(1/M, 



NJ 



2M% 



Gx{0) + — T G a {q 



AM 2 N 



^G 3 (q 2 ) + 0(1/M%) 



L N 



(45) 



^From Eq.(|45l) one can directly see that the 0(e 3 ) prediction for Ml is sensitive to sublead- 
ing order, whereas the prediction for the quadrupole multipoles E2, C2 only gives us the 
leading result. One also observes the peculiar behavior that in the case of the E2-, C2-form 
factors their q 2 -dependence is not suppressed by powers of 1/Mjv with resvect to the q 2 = 
point ! According to our knowledge this situation has not been observedtiJ before in chiral 
calculations of form factors, e.g. compare [|TT|JT~D[] E^I. Its formal origin lies in the fact that 



ii 



The notable exception has been the pseudoscalar formfactor of the nucleon which is known to 



be dominated by the light pion physics, [19|. 
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The analogous phenomenon occurs for the slope parameters of the Gi(q z 



1, 2, 3 which we 



discuss in section VII C 
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for quadrupole transitions it is not possible to write down a counter term at a lower order 
than 0(e 3 ), i.e. the same order where loops start contributing. Keeping the above caveat 
in mind that Eq. only represents the leading behavior for the E2-, C2-form factors, the 
chiral structure underlying Eq. ffidj leads us to the expectation that the E2-, C2-form factors 
only show a weak suppression at low q 2 with respect to their q 2 = point and that the 
q 2 -evolution could be very different from the corresponding (dipole behavior of the) nucleon 
Sachs form factors (^^(g 2 ), Gm (q 2 ), e.g. [|T7| , [T9|1 . On the other hand, the Ml form factor 
shows the expected 1/M suppression of its radius with respect to the photon point — quite 
analogous to the SSE calculation of the dipole form factor Gm(q 2 ) in @- This question 
will be addressed in more detail in a future communication once precision information on 
the relevant coupling constants is available. We further note that the quadrupole multipoles 
in Eq.([EJ) satisfy the long wavelength (|q| — > 0) constraint E 2 = f^C 2 - 



VII. NUMERICAL RESULTS 



To provide numerical results we need to fix the coupling constants entering the cal- 
culation. For our 0(e 3 ) calculation it is appropriate to use F n = 92.5 MeV, = 938 
MeV, 5 = 293 MeV, m n = 140 MeV, g A = 1-26, g nNA = 1-05 1H|7|. However, the 
four ^NA couplings bx, be, D\, E\ as well as the leading 7rAA coupling g\ are only poorly 
known at the moment. Until we have more information on these couplings we therefore 
proceed semi- empirically and evaluate the exact C(e 3 ) SSE results of sections II to VI with 
phenomenological input from other calculations. 

In particular, for the leading Ml coupling we use b\ = 7.7 as in Ref. f07|§, although this is 
not entirely consistent with our 0(e 3 ) calculation as this particular strength for the coupling 
only holds to 0(e 2 ) and undergoes an as yet undetermined renormalization in 0(e 3 ). The 



remaining three jN A couplings are still undetermined [^T], for now we find that we can 



parameterize them via 2 independent linear combinations C 2 , C3 and fix them from Refs. 
PJ40|1 in the next section. Finally, due to lack of more precise information at the moment, 
we use the SU(6) quark model estimate g\ = |^ for the 7rAA coupling in Eq.(|I^), which 
should be accurate to 25% (e.g. |41]]). We also note that the so called off-shell parameters 
t/i 5 6 of Eq. flT8| ) decouple from our results. 



3 The factor of 2 difference compared to ref. [37] just results from the fact that we use a different 
form for the 0(e 2 ) jNA vertex in Eq.(|l9|), in order to be consistent with the SSE formalism paper 
ref. [|l(|. The sign of bl is chosen to give a positive result for the 0(e 2 ) Ml 7./VA transition 
multipole at q 2 = 0, consistent with the previous analyses pl[| . 
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A. A best estimate for unknown couplings 



First, we note that the unknown couplings" 6 6 , D\, E\ only show up via two independent 
linear combinations in our calculation. We therefore introduce the new set of parameters 

C 2 .^l. C^^A. (46, 



We now want to utilize existing phenomenological analyses for NA transition multipoles in 
order to fix these two parameters. Most of the existing work relies on parameterizations 
of pion photoproduction in the Delta region and has focused on the determination of the 
transition multipoles at the so called K-matrix pole position, e.g. see 0. Unfortunately, 
this wealth of information is not directly transferable onto our problem because the K- 
matrix approach by construction only gives the imaginary part of the transition amplitude 
at the kinematical point where the real part vanishes. For our calculation, however, we need 
information of the transition multipoles at the T-matrix pole, where Ml, E2,C2 can be 
complex quantities even at the resonance position. This information has recently become 
available P, ^9| , ^0| . Our fitting procedure will refer to either the values of the real part of the 
transition multipoles Ml, E2 at q 2 = or to the ratio E2/M1 at q 2 = 0. We shall adopt the 
first choice for the RPI,VPI results j|0) and the second one for the Mainz results ||. The 
two parameter sets determined by fitting to ReMl(O), ReE2(0), with Re denoting the real 
part, of the RPI, VPI data will be designated A and B respectively while a further subscript 
taking the values 1 or 2 will denote the two fits used in reference [0. Specifically, using 
ifeMl(O) = 0.3GW- 1 / 2 , ReE2(0) = -0.0087Ge^- 1 / 2 we obtain 

C 2 Al = -7.80 

C 3 Al = -2.48, (47) 



while using iteAfl(O) = 0.301GW~ 1/2 , ReE2(0) = -0.048Ge1/~ 1/2 we find 

= -8.07 

C 3 Aa = -2.91. (48) 



We now consider the VPI group results. The values ReMl(0) = 0.297Gel / 1 ^ 2 , 
ReE2(0) = -0.0065GC1/- 1 / 2 give 

Cf 1 = -7.76 

= -2.72, (49) 

and finally the values ReMl(0) = OSOlGeV- 1 ^ 2 , ReE2(0) = -O.OOllGeV- 1 / 2 produce 



For the numerical analysis we utilize the scale-independent couplings D\, E\. 
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Cf 2 = -8.27 
Cf 2 = _3.3i. 



(50) 



The next step is to use the speed-plot analysi Jil of the Mainz group H . This time we need 
our expression 

EMR(0) = MW) ~ MsGM-SGM (51) 

from Eq. (|4"5|) and table |. We end up with a new parameter set, denoted by C 

EMR(0) MfflM = (-3.5 - 4.6 i)% -> C 2 C = -10.14 

-> Cf = -2.25 . (52) 



All the fits produce "natural size" values for the parameters C2, C3. Based on the above 
parameter sets we may determine an average value and a relevant error for C2, C3. We 
arrive at the values 

C 2 = -8.41 - 0.44 

C 3 = -2.73 - 0.18. (53) 



B. Predictions for the Multipoles, CMR(O) ratio 



Having fixed C 2 , C3 in the previous subsection we can now give the following predictions 
for the average values of the electromagnetic form factors and their relevant errors 

M1(0) = (308.6 - 8.9+ z28.6)Hr 3 GeV- 1/2 
£2(0) = (-6.4 - 1.7- tlQ.9)10- 3 GeV- 1/2 

(72(0) = (-10.8 - 1.7- tU.3)10- 3 GeV- 1/2 . (54) 
For the ratios we determine 



5 In view of the preliminary character of these numerical estimates we neglect the fact that the 
Mainz numbers where obtained at a pole position which would correspond to a smaller value for 
the nucleon-delta mass-splitting parameter 5. 
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ReEMR(0) = Re ) [ = (-2.52 - 0.48)%, (55) 



ImEMR(0) = Jm ^^ = (-5.24 - 0.16)%, (56) 
ReCMR(0) = Re^Mr = (-3.86 - 0.47)%, (57) 
ImCMR(0) = Im ^ffi = (-4.29 - 0.14)%. (58) 



We note that the imaginary parts of the multipole form factors (as well as the imaginary 
parts of the form factors Gi i G {1,2,3}) are the same for all fits and so we are not in 
a position to produce error bars for them. The reason is simply that the imaginary part 
of the transition amplitude stems exclusively from loop contributions which do not depend 
on C*2, C3. Of course the same is not true for the ratios of the multipoles since in this 
case real and imaginary parts of the numerator and the denominator are mixed in order to 
obtain the real and imaginary part of the ratio. We also note that our results appear much 
less sensitive to C2 compared with their dependence on C3. The reason is that this first 
parameter appears only in the G\ form factor and its contibution with respect to the other 
ones is suppressed by a power of the ratio 

In Figures §, ||, f|the g 2 — dependence of both the real and the imaginary parts of G\, G2, 
G3 is depicted, using the C, A\ and B 2 sets, respectively, for the values of the parameters 
C2, C3. As for the real part of the form factor G3 we do not want to imply that it changes 
its sign near — q 2 = 0.15GeV 2 . Given that G^{q 2 ) approaches quickly to zero even small 
corrections from higher orders can change the result for — q 2 > 0.1 GeV 2 . The resulting 
q 2 — dependence of the multipole form factors Ml, E2, C2 is shown in Figures |5], |^, |7| and 
|8] respectively. Finally the last remaining figures show the real and the imaginary parts of 
the ratios EMR, CMR. We close the present section with a remark concerning the range of 
q 2 , for which we believe our results to be meaningful. As we have already mentioned in the 
previous paragraph, the imaginary part of the amplitudes is determined entirely by a 1-loop 
diagram related to 7riV-scattering at the relatively large kinematic energy of approximately 
5 — m w ~ 150 MeV. We therefore do not necessarily expect that the C(e 3 ) calculation leads 
to a good description of the imaginary parts and point to the possibility of large higher order 
corrections in the imaginary parts of the Gi(q 2 ) form factors! This issue can only be resolved 
by a future C(e 4 ) calculation. 



C. Predictions for the Form Factors 



Though enjoying less fame than the corresponding transition multipoles of the previous 
section, the three NA transition form factors Gi(q 2 ) are actually better suited for a chiral 
analysis — the reason being that to C(e 3 ) the complete g 2 -dependence is solely governed by 
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O(e) parameters, whereas the poorly known couplings bi,b$,Di, E\ only enter in the q 2 = 
normalization of the form factors. Utilizing C 2 , C3 from the previous section one fixes 



Gi(0) = (5.50 - 0.14+ z0.20) 

G 2 (0) = (4.90 - 0.73- z7.42) 
G 3 (0) = (-1.93+ il.13). 



(59) 



We note that all three imaginary parts and the real part of Gs(0) do not depend on the 
couplings bi,b 6 ,Di, E\ and are given by the much better known leading order parameters of 
SSE. 



(60) 



We now proceed to give predictions for the slopes parameterr 1 *' 
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6- 



dt 



\t=0 5 



which can be directly calculated from Eqs. (42-44). Defining ft = ^ and N = - 397T ^2 9 F if N 
we find 
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(61) 



One can see immediately that p\ is suppressed by a power of 1/A X ~ 1/Mjy compared to 
the real photon point g 2 = 0. This behavior can also be observed in the F 2 (q 2 ) Dirac form 



factor of the nucleon (e.g. fl7jjl9fl ). An analogous result is expected for the radius of the 
Ml(g 2 ) form factor, see the discussion in section VI. 



For the slope parameters corresponding to G 2 and G 3 , t-differentiation of the J, H 
functions is needed, resulting in two divergent contributions (one non-integrable singularity 
resulting from differentiating under the integral sign and one divergent surface term from 



16 We focus on the slope parameters rather than on the transition radii in order to avoid the 
appearance of any of the poorly known couplings b±, be, Di, E± in the discussion. 
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differentiating w.r.t. the integration limit). Treating the integral as a principal- valued one 
the infinities are shown to cancel out, leaving 



p\ 



M 



N 



2p - -7T/i 2 
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3 P 
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1 + 8fi 2 



-3/i 2 log(- 



p\ = -2irN 



M, 
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16 , , 15 9x fir , A 1 
^ + y/^ 3 - (1 + yM 2 ) U " * ° g( ~ + ~ 1 




It is important to note that this is a quite peculiar result. One can see immediately that these 



two slope parametersG-3 scale as A ° 



i.e. they show no suppression by a heavy scale 



with respect to the q 2 = point! This behavior is completely analogous to the one expected 
for E2(q 2 ), C2(q 2 ) (cf. section |VT]), but has not yet been observed in chiral calculations of 
baryon form factors. The first moment with respect to q 2 of these "quadrupole-like" form 
factors G2(q 2 ), G%(q 2 ) is entirely given by the small scales m^, 5 and leading order coupling 
constants — allowing, in principle, a large variation of the form factors at low q 2 compared 
to their q 2 = point. However, we want to emphasize that Eqs. (JH],^) only represent 
the leading results for the three slope parameters, more cannot be learned from a 0(e 3 ) 
calculation. At present, we have no way of knowing the size of the 0(e 4 ) corrections, but in 
principle it is possible to calculate them at a later staged. 

Finally, we can obtain the chiral limit (fi — > 0) behavior of the slope parameters from 



Eqs.flU, |D 



3 , 

h log 2 

2 3 & 



log A* 



+ 0(ai) 



pI = 4"*^ + °^ 



pI 



-2nAf 



M 



N 



7T 



2 +^0og2 
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+ 0(p) 



(63) 



Thus, the real part of p\ diverges logarithmically in the chiral limit while its imaginary part 
remains finite. On the other hand, the real part of p\ vanishes altogether in the chiral limit, 
while its imaginary part remains finite. As for p|, it is the real part that remains finite in 
the chiral limit while its imaginary part diverges logarithmically. The numerical values for 
the slope parameters are also collected in Table [II] and since they take contributions only 
from the loops are independent from the values of C 2 , C 3 . 



17 The same of course holds for the corresponding radii. 

18 At first one would need a precise determination of b\ and the subleading corrections to the -kNA 
vertex, as these couplings would enter the new loop diagrams of 0(e 4 ). 
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VIII. SUMMARY AND OUTLOOK 



The main features of our investigation can be summarized as follows: 

• The A — > N'-f* amplitude has been calculated to 0(e 3 ) in the "small scale expan- 
sion" formalism. The pertinent counterterm contributions, relevant 1/M corrections 
to lower order coupling and all loop diagrams allowed to this order have been analyzed 
systematically. It was found that the g 2 -dependence of the transition is generated by 
two tyN— , Ti A— loop diagrams and that two counterterms are needed to renormalize 
the result. 

• We have obtained an estimate for presently poorly known couplings from existing 
phenomenological analyses and then given predictions for the complex transition mul- 
tipoles Ml(g 2 ), E2(q 2 ), C2(q 2 ) as well as the complex multipole ratios CMR(O), 
EMR(g 2 ), CMR(g 2 ). While matching the 0(e 3 ) results to the corresponding expres- 
sions of the multipoles as functions of q 2 we found that the radii of the quadrupole 
form factors are not suppressed by a large scale compared to the q 2 = point, indi- 
cating a possibly rapid variation at low q 2 . Further study of this interesting result will 
continue once more precise information on the necessary couplings is available. 

• Three appropriately defined electromagnetic NA transition form factors Gi(q 2 ), i = 
1,2,3 have been identified from the 0(e 3 ) A — > N'-y* calculation. The longitudinal 
form factor Gs(q 2 ) had to be rescaled compared to existing definitions in the literature 
in order to achieve consistency with the power-counting of SSE. Gi(q 2 ) could be found 
to subleading order, whereas for G 2 (q 2 ),G3(q 2 ) we obtained the leading result. 

• The three Gi(q 2 ) form factors are found to be complex quantities due to the 7riV-loop 
diagram (i) in Fig.l. Their entire ^dependence is controlled by relatively well-known 
0(e) parameters — making the Gi(q 2 ) form factors the prefered testing ground of chiral 
symmetry in the radiative iVA-transition. The corresponding slope parameters have 
been calculated and their chiral limit behavior has been discussed. 

The numerical results of this work have to be considered preliminary due to the presently 
poor state of knowledge on coupling constants in the SSE formalism. However, this situation 
is going to improve in the near future and we will at a later point revisit our 0(e 3 ) analytical 
results and try to improve upon the numerical accuracy of our predictions. 
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APPENDIX A: INTEGRALS 



The A 1 ^ (fl, M 2 ) functions used in Eq.(|36[) are denned as Euclidean integrals in dimen- 
sional regularization 

A^^M?) = / oP da / f^r— 5 ,\ E > 7777- ■ (Al) 

v ' * ; Jo J (2vr) a! (-1% - \a 2 + aQ - M?)* V ; 

To calculate them, first of all, we express the denominator using the Schwinger proper time 
representation 

f° dxx n - l e- Ax = ^ A>0 (A2) 
Jo A n 

with r(n) the well-known F function. The constraint A > is explicitly fulfilled for the 
Q < case while for Q > we analytically continue our expressions. Finally using the 
tabulated integrals 

jf x^e-^dx = (2/3)- W2) r» e*p(^)£U (j^yji) (A3) 

r°° rt ,„/, vi,, 2 1 - /3 -fr(/3) . A i /i/^ y + ^ + i 

jf e-«r^D M W = r( + p^T) (» + ( 5 ■ f i — g~ ; 7Tl) (A4) 

with D_ v (x) the parabolic cylinder function and F(a,b;c;x) the hypergeometric function, we 
find explicitly 

_fW _ Ag\ _L_ (-y ^-Q 2 arccos (--fe) 

d " + 32vr4 1 %^ + 167r 2 x | +v ^r^i n ^ 1 | + ^TT) 



A 301 1 n ( M?\ 1 [ - ^M 2 -^ 2 arccos (--|) 

d(d-l)~ 3^+ 288^ 31n ^J + 48vr 2 X | + ^3^1n(--| + ^ 



i / )i iri \ o 7 I arccos 

^.4-™ = ?L + -J^ (7 + 3 in 4) - A x ( 



1 I / I arccos ( — 4r 

A 300 = L_ I V^F^ V M< 

16tt 2 I i i n ( I ,/» 2 7 



d3 A 300 = 1^300 ) (A5) 
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where the factor L carries the infinity in dimensional regularization and is defined in Eq. 
(j39[). It is clear that A 300 is finite while A 310 and A 301 diverge. The integrals in Eq . (|37|) are 
split between [0,xq] (trigonometric branch) and (xq, 1] (log branch) where xq is the large 



root of Q 2 = M 2 , 



-t + sjt 2 + 4m2|q| 
x o = iTi - !^ 



2|q| 

with t = 5 2 — |q| 2 to the order we work. In the calculation the functions are defined with 
Q = 5x for the Ntt intermediate state loop diagram and Q = —5x for the An one. This 
entails (cf. Eq.(X5)) that the logarithms have negative arguments (and therefore there is 



an absorptive piece of the amplitude) for the Ntt diagram alone, as expected from general 
considerations (s = M\ > (Mi + m n ) 2 for i = N only). The imaginary parts are computed 
analytically. 

For completeness, we show the correspondence between our A functions and the J func- 
tions appearing in the literature (e.g., 



a 200 (q,m 2 ) = J (n,M 2 ) 
A 20 \n,M 2 ) = dj 2 (n,M 2 ) 

A 3 ^M 2 ) = \-^-U^M 2 ) 



A 30 \n,M 2 ) = ^j 2 (n,M 2 

Z M 2 



d d 

_/./<> i / 

t 



A™(Q, Mt) = -~J2(P, Mt) = Ji(fl, Mt) . (A6) 
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TABLES 



TABLE I. The ANj form factors: a = - 97lN ^ N and [I] is short hand notation for 
[9aIn + X9iIa] , with / G {F, J, H} 



Source 


Gi 


G 2 


G 3 


Loops 




4a[J] = 52(g 2 ) 


4aS[u] = Uq 2 ) 


A {2) 


bi 
2 


61 










-61 - 26 6 










n 


n 
u 




4M JV ' D 1 







A { V 










TABLE II. The slope parameters [fm 2 ] 


of the form factors Gi(q 2 ) 


1 


2 
Pi 


pi 




2 

Pi 


1.94 + 0.55 i 


10.50-6.16 i 




-4.12 - 4.48 i 
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FIGURES 




(i) 0) 



FIG. 1. Loop diagrams: single(double) solid lines denote Nucleons(Deltas), respectively. 
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FIG. 6. The real part of the form factor E2 with C2, C3 given by the sets A\, B2, C. 



-0.005 



-0.01 



-0.015 if 



-0.02 



> -0.025 
O 

-0.03 



-0.035 



-0.04 



-0.045 



~i r 



1 

ReC2_A1 
ReC2_B2 + 
ReC2 C □ 





□ o 



o 

□ o 

□ 

□ o 
□ o 



□ 



□ 



□ " 

□ o 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

-q"2[GeV"2] 

FIG. 7. The real part of the form factor C2 with C2, C3 given by the sets Ai, B2, C. 



34 



-0.01 



-0.012 



^ -0.014 



CD 



-0.016 



-0.018 



-0.02 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

-q A 2[GeV"2] 



0.2 



FIG. 8. The imaginary part of the form factors E2, C2 with C 2 , C3 given by the sets Ai, B 2 , C. 



-0.5 
-1 

-1.5 
-2 

-2.5 



-3.5 [*■ o 

T □ o 



-4.5 
-5 
-5.5 



~i 1 

ReEMRAI o 
ReEMR_B2 + 
ReEMR C □ 



□ □ 



□ □ 



□ □ 



□ □ 



□ □ 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

-q"2[GeV"2] 

FIG. 9. The real part of the EMR ratio with C 2 , C 3 given by the sets Ai, B 2 , C. 
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FIG. 10. The imaginary part of the EMR ratio with C2, C3 given by the sets A\, B2, C. 



-3 



-5 



~i r 



~i 1 1 r 



~i 1 

ReCMRAI o 
ReCMR_B2 + 
ReCMR C □ 



o n 
o a 
o □ 



° ° □ 



5 o> 



_1 L_ 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

-q"2[GeV"2] 
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FIG. 12. The imaginary part of the CMR ratio with C2, C3 given by the sets Ai, B2, C. 
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